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Abstract 

We study the supersymmetric Wilson loops in the four- dimensional = 4 super Yang- 
Mills theory in the context of AdS/CFT correspondence. In the gauge theory side, it 
is known that the expectation value of the Wilson loops of circular shape with winding 
number fc, Wk{C), is calculable by using a Gaussian matrix model. In the gravity side, 
the expectation value of the loop is conjectured to be given by the classical value of the 
action ^ds for a probe D3-brane with k electric fluxes as {Wk{C)) = e"'^^^ . Given such 
correspondence, we pursue the interpretation of the matrix model eigenvalue density, or 
more precisely the resolvent, from the viewpoint of the probe D3-brane. We see that the 
position of an eigenvalue appears as an integrated flux on the D3-brane. In the course 
of our analysis, we also clarify the boundary condition on the D3-brane in terms of the 
Wilson loop. 
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1 Introduction 

Duality between string theory and gauge theory has been an important concept for the- 
oretical particle physics. In particular Wilson loop would play a unique role there, since 
originally it has been introduced in the context of stringy behavior of the strong coupling 
gauge theory, i.e., the area law behavior in the confining gauge theory. 

In the recent development of the duality, the "area law" of the Wilson loop has been 
rediscovered [1] in the context of the AdS/CFT correspondence. Hence in AdS/CFT, 
Wilson loop is useful for concrete realization of the duality. There, we usually begin with 
the N D3-branes located on top of each other and put another D3-brane in parallel with 
them at a large distance. Then we evaluate the amplitude for the propagation of a string 
stretched between the N D3-branes and the single separated brane from two different 
points of view, from gauge theory and from gravity picture. 

From the viewpoint of the gauge theory, the stretched string corresponds to a bi- 
fundamental matter of the SU(iV)xU(l) gauge theory. By evaluating the amplitude for 
the string propagation along a loop C on the isolated brane, we can introduce a Wilson 
loop for C in the bi-fundamental representation of the SU(iV)xU(l) gauge group. 

In order to discuss the gravity side, we replace the N D3-branes with the near horizon 
geometry of the extremal black 3-brane solution. Then the isolated brane is recognized 
as a single probe D3-brane in the near horizon region and we evaluate the propagation 
amplitude for the string attached to the loop C on the probe D3-brane. 

The AdS/CFT correspondence claims that these two different points of view give the 
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same result. Therefore we reach the following conjecture: 

—trPexp ( [ dsiiA^x^" + \x\9i<i>i))\ = I rfXe"^^*'-^ . (1.1) 

N 'I CFT JC 

Here on the left hand side, not only the gauge field but the scalar fields are included since 
a string attached to D-branes is coupled to all of them. Also, the U(l) part of the Wilson 
loop has been omitted since we discuss the limit in which the U(l) brane is separated by a 
large distance from the remaining D3-branes and thus its dynamics decouples from the 
SU(iV) part we are interested in. Actually in the gravity side we consider the string world 
sheet which is attached to the loop C on the AdS boundary. Hence the right hand side 
only takes account of the SU(A^) contribution in terms of the gauge theory side. Also we 
need to take care of the definition of the functional S'string on the right hand side. In the 
paper [2] it has been proposed that in addition to the usual Nambu-Goto type action, we 
need to add appropriate boundary terms which eliminate divergence due to the infinite 
scale factor from the vicinity of the AdS boundary. Introducing these boundary terms was 
recognized as performing the Legendre transformation to change the boundary condition 
of the world sheet. 

The conjecture (11. ip has been checked mainly in the case of the straight Wilson line or 
the circular Wilson loop which preserves some global supersymmetryj^ In the study of the 
circular Wilson loop, a Gaussian matrix model plays an important role. It was proposed as 
a technical tool for summing up all the planar ladder diagrams [4] and a further argument 
for the matrix model including the non-planar diagrams was given in [SjH In the more 
recent paper [7], an argument based on the topological twist and localization technique 
is discussed. A relation to the Gaussian matrix model has also been discussed by using 
the mirror symmetry and the geometric transition [8] . 

In the study of the validity of the conjecture (11.11) . an important work has been done in 
the paper [9]. There the authors considered a circular Wilson loop with winding number k. 
This winding number corresponds to the string charge because k stretched strings between 

D3-branes and the single D-brane mentioned above yield the Wilson loop with winding 
number k. In the gauge theory side the expectation value of the operator can be evaluated 
by the matrix model of the same Gaussian action as in the single winding case. On the 
other hand, in the gravity side, they considered a spike D3-brane solution carrying a non- 
trivial electric flux by k units of the electric charge, because the electric charge on the 
D-brane corresponds to the string charge. In the gravity side the computation is essentially 
the same as the right hand side of (11.11) with the string action being replaced by the D3- 
brane actionjf] In [9], it was found that the relevant D3-brane solution has the geometry 



^ For a recent quantitative test of a similar conjecture to in the case with less symmetry, see the 
paper [3]. There a similar relation in the case of finite temperature DO-brane system is tested by using 
the Monte Carlo simulation. 

^See also [6] for a discussion on the relation between the Wilson loop and the matrix model in the 
beta-deformed super Yang-Mills theory. 

^ It has been proposed that the D3-brane actually corresponds to the Wilson loop in the symmetric 
representation [10], and has been argued that in the strong coupling limit of the gauge theory, a multiply 
wound Wilson loop and a symmetric one give an identical expectation value. See the paper [11], for these 
issues. We however consider that there is a subtlety here. See the concluding section. 
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of AdS2xS^ where the radius of the is given by the parameter k, = k^/X/{AN). This 
suggests that the large- iV hmit with fixed k in the gauge theory side is a quite interesting 
hmit. Actually it was found that in this limit, the expectation value of the Wilson loop 
computed by using the Gaussian matrix model agrees with the prediction of the gravity 
side. A very important point in this work was that the large-iV limit is not the planar 
limit but the result contains a class of the non-planar contributions since the parameter 
K depends on N inversely. 

One of the interesting developments which followed the paper [9] is the study of the 
geometric aspects of the eigenvalue in the Gaussian matrix model. Stimulated by the 
case of a local operator [12], the correspondence between the eigenvalue distribution and 
a certain aspects of the geometry has been discussed in [13,14]. In general, gravitational 
interpretation of eigenvalues in matrix models would be an interesting and essential prob- 
lem for gauge/string duality, or even for nonperturbative formulation of string theory. For 
example, in the IIB matrix model [15] the eigenvalues of matrices are interpreted as the 
space-time points and their dynamics is discussed as emergent geometry [16]. Therefore, 
even if we concentrate on a particular eigenvalue of a specific matrix model, it would be 
intriguing to find its gravitational interpretation clearly in the context of the AdS/CFT 
correspondence. 

Information on the eigenvalue of a matrix model is well packed in the expectation 
value of an operator called the resolvent. In this paper we thus consider the resolvent of 
the matrix model from the viewpoint of the probe D3-brane. Our main goal is to discuss a 
gravitational interpretation of the eigenvalue distribution, or more precisely the resolvent, 
in the Gaussian matrix model in the presence of the Wilson loop with a large winding 
number k. In particular we are interested in results which follow directly from the basic 
correspondence (11. ip using the D3-brane action. This is in contrast to the "bubbling 
geometry" approach taken in [12-14]. On the way of the analysis, we also identify what 
kind of boundary conditions should be imposed on the D3-brane configuration so that it 
will correspond to the Wilson loop with winding number k. 

The rest of the paper is organized as follows. In section [21 we review the matrix model 
computation and discuss the resolvent in the 1/A^-expansion with fixed k. From this 
result, we will see what kind of quantity in the gravity side should reproduce the matrix 
model resolvent. In section [3l we begin with a brief review of the setup and the results 
of [9] . Next we perform the computation in the gravity side which yields the same result as 
in the matrix model, and then we find a gravitational interpretation of an eigenvalue. In 
the course of the analysis we will give some interesting aspects of the work [9] , especially 
identification of appropriate boundary conditions for D3-brane configurations. We then 
conclude the paper in section H] with some discussions. A couple of appendices are devoted 
to filling in some technical details in the main part of the paper. 

2 Multiply-wound Wilson loops in gauge theory 

In this section we discuss the gauge theory side and its relationship to the matrix model. 
In subsection 12.11 we start with a brief review of the expectation value of the Wilson 



3 



loop in the large-iV matrix model. We also mention the eigenvalue distribution derived 
in [13] by solving the saddle point equation of the matrix model. For our present purpose, 
however, it will turn out to be more useful to discuss the resolvent through the Laplace 
transformation. In subsection 12.21 we discuss the Laplace transformation that motivates 
us to proceed to the gravitational interpretation. 

2.1 Matrix model calculation 

We consider the Wilson loop in the four- dimensional J\f = 4 super Yang-Mills theory, 
which is defined by 

Wk{C) = ^trPexp j ds(iA^{xis))x^'{s) + . (2.1) 

Here A^{x) {fi = 1 . . . 4) and (i = 1 . . . 6) are the gauge field and the scalar fields, 

respectively. The coordinates span the Euclidean four-dimensional space and 6^ is the 
coordinate on the unit . We consider the operator in which the path C : {x'^(s) | < 
s < 2A;7r} is a multiply winding circle and 6^ is constant. The subscript k on the left hand 
side indicates that the loop C goes around the circle {x^{s) | < s < 2ti} k times. 

The expectation value of the operator fl2.ll) is calculable by means of the Gaussian 
matrix model [4, 5, 7] : 

e'^L^^U-k")^f{k',N), (2.2) 

(2.3) 

Here k' = ^/X/{4N)k and Lt\0 = (e«C""/^!)c^"/c^C(e~^C"^°) is the Laguerre polyno- 
mial. We introduced the function f{k',N) for later convenience. Hereafter, the expecta- 
tion value of the Wilson loop (12.11) is always understood as this matrix model expectation 
value. 

In [9], the large- limit of f{k', N) with fixed k = ky/X/{AN) was derived by using a 
differential equation. Here we follow their steps and introduce J-'{k, N) = —N~^ log f{k', N) . 
Then JF satisfies the following differential equation: 

{d^^y - + 39,^) - 16(1 + /t^) = . (2.4) 

kJS 

In the large- iV limit with fixed k, the differential equation (12. 4p can be solved perturba- 
tively and the solution is given bjfl 

^ The constant term can be fixed, for example, by comparing the small k, limit of (|2.5p with the 
modified Bessel function. For the relation between the function /(fc', N) with small k, and the modified 
Bessel function, see the discussion around (|2.28p . 



:^^(^)>Ar.4SYM=(^tre^*^ 



matrix model 



Z j N 

Z^ I dMe-"-^'^^'' 
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= ±2 (^kVTTk^ + arcsinhfi;) + ^ log «Vl + + log(327rA^3)^ + 0{N~^) . (2.5) 

We take JF_ since it dominates in the large- iV limit. 

Now let us see that the leading term in (12. 5p originates from an eigenvalue apart from 
the cut of the standard semi-circle distribution of the Gaussian matrix model [19]. In 
terms of eigenvalues, (12. 2p can be written as 

{Wk{C)) = ^J Hdrmexpi-Ves), 



N-l N-1 N-1 



Vef[ = N V{mi) — log(mj — rrijY + NV{mj\f) — log(m7v — mjY — kmi 

i=l i>j j=l 

V{m) = ^m^, (2.6) 
A 

which tells us that the Wilson loop gives rise to force only on a single eigenvalue (say 
itin)- Therefore let us introduce the eigenvalue distribution as 

p(m) = p''^\m) + —p''^\m), 
p''°-*(m) = — (5(m — mj), p^^\m) = 5{m — m]\j). (2.7) 

i=l 

Then can be rewritten as 

dmp^^\m)V{'m) — / dmdm p^^\'m)p^^\m')log{m — m'Y 

+ N J dmp''^\m)V{m) — N J dmdm! p''^\m) p^^\ra') \og{ni — m'Y — k J dmp''^\m)m. 

(2.8) 

By considering the variation of p^^\m), the saddle point equation reads 

V■^,„^ _ 2 Un'^^ - 1 Ln-f^ = 0. (2,9) 

J m — m' N J m — m' 

In order for the subleading distribution p^^\m) to make sense, we need to discuss the 
distribution p^^\m) to the same order. Hence we further expand it with respect to 1/N 
as 

pio)(m) = p(°'°)(m) + lp(0'i)(m) + C(Ar-2). (2.IO) 

Here we note that the distribution functions p^^'^\m) and p^^'^\m) satisfy the following 
conditions: 

c/mp(°'°)(m) = 1, I dmp^°'^\m) = -1, (2.11) 



which can be seen from (12.71) and the expansion fl2.10p . Hence, formally, p^^'^\m) corre- 
sponds to the distribution of N eigenvalues while {1/N)p^^'^\m) subtracts a single eigen- 
value. By taking terms of order in fl2.9l) . we find that p^'^'^\m) satisfies the saddle 



point equation for the semi-circle distribution 27rp*^°'°^(m) = y y ~ ^rn'^ with support 

Next, before discussing p(°'^)(m), we shall study the saddle point equation which will 
be derived by considering the variation of p^^\m). In order to solve this, it is reasonable 
to assume that mA? is isolated from the other — 1 eigenvalues due to the extra force 
coming from the Wilson loop and to make an ansatz p^^\m) = 5{m — m^) with outside 
the range (— \/A, VX). Then the saddle point equation for positive m* becomes 



y'K) - \^V'{m,) - ^V'{m,y _ - A = 0, (2.12) 
from which we find 

m, = yXVl + (2.13) 

which is indeed outside the cut|f| 

Finally, we discuss p^^'^\m). Let us first study the behavior of p^'^'^\m) in the limit 
/t — s> oo by a physical argument without using the saddle point equation. In this limit, 
the isolated eigenvalue goes to infinity and its effect on the remaining — 1 eigenvalues 
should vanish. Then the resulting distribution function p^^\m) becomes the semi-circle 
distribution for — 1 eigenvalues which is given by 

/V — 1 2 1 

lim p(°)(m) = — — X VA^_i - m2 + 0(Ar-2) ^2.14) 

K^co N IT AaT-I 

= v/A^ - ^l^=L= + 0{N-'). (2.15) 

Here we have introduced Atv-i = {N — l)/N x A and the overall factor (A^ — 1)/A^ on the 
right hand side of (12.141) reflects the fact that p^^\m) is defined with the overall factor 
1/N instead of 1/(A^ — 1) (see (12. 7p ). So, from this simple argument we find that, in the 
limit K CO, (1/N)p^^'^\m) is given by the second term of (12.151) (including the minus 
sign) . 

For a generic value of k, we need to solve the saddle point equation for p(°'^)(m), i.e., 
the terms of 0{N~^) in (12.91) . which is now given by 

^(0,1) f^/^ 1 

dm'P ^ = 0. (2.16) 

m — m' m — m* 

Physically, this equation shows that the subleading distribution p^^'^\m) is determined 
by the following manner: the repulsive force from the isolated eigenvalue at m = 



''This is the unique solution for fc > 0. For negative k, we would have the same solution with the 
opposite sign. In this paper we assume that fc > without loss of generality. 
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changes the subleading distribution p^^'^\m) of the — 1 eigenvalues from its form in 
the hmit k, ^ oo discussed above. This change gives rise to additional force among 
the — 1 eigenvalues themselves. Then the distribution p^^'^\m) is determined by the 
force balance condition among these effects. Therefore, p^°'^'*(m) describes both effects 
of subtracting a single eigenvalue and distortion on the distribution of the remaining 
eigenvalues due to the isolated eigenvalue. It is evident that they are of 0(1/N). Since 
p(°'^)(m) subtracts an eigenvalue, it is natural to assume that it has its support only inside 
the range or at m = m^. This is because the total eigenvalue density p{m) 

must be positive everywhere. In fact, as we have shown above, the support of p^'^'^\m) 
in the limit k ^ oo is given by (— V^, V^)- Furthermore, since the distortion due to the 
isolated eigenvalue should reach the whole distribution of the remaining eigenvalues, it is 
natural again that the support of p^^'^\m) is {-y/X,VX) for a generic value of k. In the 
rest of this subsection we assume this propertyj^ 

Let us solve the saddle point equation fl2.16p based on this physical argument. For 
this purpose, we introduce the resolvent R^^'^\m) as 

Ri^'^Hm) = - fdm'P-!-:^. (2.17) 
^ ' N J m-m' ^ ' 

Then (m) is determined by the following conditions: 

• The real part of i?^°'^)(m) at the support of p^°'^)(m) is given by (12.161) . 

• For the generic value of k, the resolvent has one cut (— -\/A, \/~\) while it has no pole. 
This is due to the assumption for the support of p'^'^'^\'m) mentioned above. 

• In the large-m limit, the condition _R'^°'^)(m) — — j^^ should be satisfied. This is 
due to the second equation of (12.111) . 

By taking account of these conditions, we can fix the resolvent and the eigenvalue distri- 
bution completely as 

N\ m-m^ {m - m^)\/im? - \ J 
p(0'i)(m) = -;^fi?(°'^)(m + ze) -i?(°'^)(m-«e)') (2.19) 

' -v^<m< v^) («^0), 



TT (m - ■m^)\/\ - m? ' (2.20) 

-5{m-^f\) (/t = 0). 

Here in the case of k = 0, we have the delta function distribution which exactly cancels 
the another subleading distribution p^^\m). This is indeed the expected behavior since 



^ In the next subsection and appendix lA.ll we give a rigorous discussion of this model by using the 
orthogonal polynomial method without making any assumption. In particular, in appendix lA.il we derive 
the explicit form of the resolvent independently, which also provides a proof of this property. 
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for K = 0, the distribution should reduce to the leading one p(°'°)(m). It is also easy to 
check that in the limit k —>■ oo, the above distribution shows the behavior of (A^ times) 
the second term in fl2.15p . 

Having studied the eigenvalue distribution, we now consider the saddle point value of 
Ves- The leading term of order N"^ is just canceled by the same term from the denominator 
Z. The sub leading term of order which comes from the first line of (12. 8p is given by 

^a),ist line _ /" ^^p(o,i)(^) /^y(^) _ f dm' p^'^'^^^m') log(m - m')A . (2.21) 



By using the saddle point equation, it is easy to show that the combination in the round 
bracket, which is obviously /c-independent, is also m-independent. From this, we find that 
(I2.2ip is just a fc-independent constant, since m-integral of p^'^'^\m) is /c- independent. 
Next, the order N term coming from the second line of (12. Sp is given by 

^a),2nd line ^ ^ (^-2kVI + - 2arcsinhK + constant) , (2.22) 

where the last term is a ^-independent constant. We can check that the constant term in 
(I2.22P is canceled by (12.211) due to the relation p'^'^'^^m) + p^^\m) = for k = 0, which 
has been mentioned after (I2.20p . Therefore, we have derived the first term of (12.51) with 
confirming again that we should take the minus sign there. An interesting point in this 
analysis is that, since we consider the large k ^ N limit, the effects of the Wilson loop 
on the saddle point are not negligible; the standard semi-circle eigenvalue distribution is 
corrected by the 1/A^ terms and these corrections affect the leading behavior of the Wilson 
loop expectation value. In particular, the important contribution of (I2.22p originates from 
the isolated eigenvalue 

2.2 Resolvent through the Laplace transformation 

The eigenvalue distribution discussed in the previous subsection is a kind of master field 
for the correlators in the presence of the Wilson loop. If AdS/CFT can be regarded 
as realizing the idea of the master field, it is natural to expect that there should be 
a counterpart of the eigenvalue distribution in the gravity side. Here we note that in 
general multi-matrix models like the four-dimensional A/" = 4 super Yang-Mills theory 
would not have a master field, but in a certain BPS sector like circular Wilson loops a 
kind of master field may exist due to higher supersymmetries by which the system can 
be effectively described by a one-matrix model. As we have already used in the previous 
section, the eigenvalue distribution is given by the imaginary part of the resolvent, R{z) , 
which is more tractable as a complex function, and the position of the eigenvalue is the 
pole of R{z). We thus study how the isolated pole arises when we discuss the resolvent 
based on its standard definition: 



^Such isolated eigenvalues play a crucial role in nonperturbative effects in noncritical string theories 
[17, 18]. It would be interesting if the Wilson loop with winding number of 0{N) is shown to play a 
similar role in critical string theories. 
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= dpe-P'^ J rfM^^tre^'^e-^W , (2.24) 
Zv= [ dMe-^^^^\ (2.25) 



Here V{M) is a potential whose saddle point we are interested in. Such a question may 
be interesting even purely from the matrix model viewpoint. What is more, it will guide 
us to a gravitational interpretation, as we shall see. 

In the integral representation (12.241) . we should pay attention to the range of the 
variable p . Since we are interested in the saddle point which is defined solely by the 
potential V{M) , we need to drop effects of the inserted operator tre^^ on the saddle 
point. This means that although the upper limit of the p- integral is set to be infinity, 
we still assume the relation p <^ N . In order to make this point clear, let us recall the 
case of the Gaussian potential, V{M) = {2N/ A)tr . In this case, the resolvent (12.241) 
is given by 

i?(^)Gau..ian = dpe'^^ f {p' , N) , (^p' = ^J^p^ . (2.26) 

The function f{p', N) , which is defined in (12.21) . satisfies the following differential equation: 

p'dl,f{p', N) + 3dp>f{p', N) - p'{m + p'')f{p', N)=0. (2.27) 

By using the relation p' = a/ A/ {AN)p and by taking the large- limit with p ^ N , (12.271) 
is reduced to the Bessel differential equation and then we have 

/(p',iV)-^-^/i(pv^), (iV^oo, p«iV). (2.28) 

Here Ji is the modified Bessel function and the overall constant can be fixed by requiring 
f{p', N) — 1 in the limit p —>■ 0. Inserting this into (12.261) . and performing the p-integral, 
we find the resolvent for the Gaussian potential: 

rfpe-P^— ^/i(pv^) = jiz - ^/^^) . (2.29) 

Next we turn to the case with the insertion of a Wilson loop with a large winding 
number k of 0{N). In this case, we need to take account of effects of the Wilson loop on 
the saddle point, i.e., we should consider the following potential: 

g-y(M)^l^^gfcMg-^trM^_ (2.30) 

Hence the resolvent in the presence of the Wilson loop with winding number k, which we 
call Rk{z), is defined by 

/•oo 

Rkiz) = / dpe-^'W{p)k, (2.31) 
Jo 
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W{p), ^j-J ciMltre^^^ltre'^^e-^*^^^ , (2.32) 

Z, = j dM^^tre'^^e-"-^''^'" . (2.33) 

Here we note that the computation in ( I2.3ip and (12.321) should be done in the following 
order: first we calculate the "two loop correlator" (I2.32p for finite and then we take 
the large-iV and large-fc limit with k and p kept finite. Finally we perform the Laplace 
transformation (12.311) . 

An explicit computation of the two loop correlator W{p)k is given in appendix lAl by 
going to the eigenvalue integral and using the orthogonal polynomial method. In the 
eigenvalue integral, we observe that the following decomposition, 

N N 

tr tr e*=^^ = J2Y1 e^"''e'="^ + ^ 6^^'+'=)"^' , (2.34) 

i=l jj^i i=l 

naturally arises and then these two terms are analyzed separately. A physical interpreta- 
tion of these terms will be discussed shortly. Here we just call the expectation value of 
the first term w{k,p), and the second w{k + p): 

W{p)k = w{k,p) + w{k + p). (2.35) 

The explicit results of w{k,p) and w{k + p) are given in (1A.30I) and (1A.31I) : 

wik,p) = fip,N)--e2 ^ .(1) ' (2.36) 

ij ^N-l\ ^ ) 

1 f{k'+p',N) 

"(^ + ^^ = iV fik',N) • ^'-'^^ 

The first term on the right hand side of (I2.36P gives the resolvent which is identical to 
the one for the Gaussian potential (12.291) in the large- A^ limit after the Laplace transfor- 
mation. The second term gives the 0{N~^) resolvent (12.181) . as shown in appendix lA.li 
Here note that in the previous subsection, we assumed that the 0{N~^) resolvent (12.180 
has the cut at (— V^, V^) and has no pole (see the second condition after (12.170 ). How- 
ever in appendix lA. II we explicitly derive the resolvent (I2.18P by taking the large- A^, fixed 
K and p limit for the second term in (12.361) and performing the Laplace transformation, 
which justifies the assumption. 

Next let us consider the Laplace transformation of (I2.37P . We should take the large- A^ 
limit with assuming the relation p -C A^ before the Laplace transformation. Furthermore, 
this term depends on the additional parameter k and the parameter region is chosen so 
that the combination k = ky/X/ (AN) will be kept finite. The relevant terms in this limit 
are already given in (12. 5p . We just define k = {k+p)y/X/ (4A^) and then (I2.37P is evaluated 
in the large- A^ limit with k fixed and p -C A^ as 



1 f(k' + p',N) 1 

exp 



N f{k',N) N 



N{J^4R,N) -J^^{k,N))^ (2.38) 
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1 




(2.39) 



N 



exp 



By performing the Laplace transformation, we obtain the following contribution to the 
resolvent of 0{1/N): 



This is exactly the pole corresponding to the isolated eigenvalue f l2.13p found in the 
papers [13,19]. Since the other 0{N~^) contribution, the second term in w{k,p), is 
shown to give only the cut, these two 0{N~^) contributions therefore give clearly different 
contributions to the resolvent or eigenvalue distribution. 

We have thus succeeded in identifying which correlator gives rise to the isolated pole, 
but, of course, the result fl2.40p itself is not quite new since we have just calculated the 
position of the isolated eigenvalue by a different method. However, the above procedure 
implies the following interpretation of the Laplace transformation and the resolvent. 

First note that, as we shall review in subsection 13. 11 the Wilson loop tr e^^ corresponds 
to the D3-brane with k units of the electric charge in the dual gravity picture in the 
AdS/CFT correspondence (11.11) . Then the two loop correlator (I2.32p may be regarded 
as treating the D3-brane with string charge k as "a part of the background" and then 
introducing additional small number of string charge p , by which the total background is 
probed. 

Now we may advocate the following physical interpretation of each term in the de- 
composition (I2.34|M the second term can be regarded as a situation where the additional 
string charge p is dissolved in the D3-brane with charge k, now altering the number of 
the charge from k to k +p. On the other hand, the first term may correspond to the case 
in which the additional string charge p is described by a probe D3-brane with charge p 
that is separated from the "background" D3-brane with charge k, but has an interaction 
with it. At the leading order of the large- iV limit, we may expect that the interaction 
between them is not taken into account and that the first term of (12.361) describes two 
non-communicating D3-branes. Actually the first term in (12.360 is blind to A; as a result of 
cancellation of the effect from the D3-brane with charge k between the numerator and the 
denominator in (I2.32p . while the second term in (12.360 describes the interaction between 
these D-branes in higher order of the string coupling constant. It then seems natural 
that the first term of (I2.35P reproduces the semi-circle distribution, which would carry 
information on the AdSsxS^, while the second term gives the pole, which is information 
of the D3-brane with k charge, in the leading contribution with respect to 1/A^ of each 
term. 

Based on the above derivation of z^, , we may also expect that, in the gravity side, 
z^: appears as an object which is conjugate to the string charge in some sense. We will 
further investigate these points from the gravity side in the following section. 

^In [19], a similar decomposition of the exponential operator is discussed in terms of world sheet 
nonperturbative effects, in a slightly different context. 




(2.40) 
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3 Matrix model resolvent from D3-brane picture 

In this section we pursue the gravitational interpretation of the matrix model resolvent. 
In particular we concentrate on the pole corresponding to the D3-brane with string charge 
k . We will discuss the cut in the concluding section. 

Since we found that the isolated pole is derived by the Laplace transformation of the 
expectation value of the operator tr e^'^+P^^'^, it is rather easy to derive the pole itself in the 
D3-brane picture once the correspondence between the operator tr e^'^^ and the D3-brane 
solution with k flux is established. Indeed the expectation value of the operator tr e^'^^^^'^ 
in the large- and large-fc limit, or equivalently the function J^{k, N) is exactly reproduced 
by using the D3-brane solution as reviewed in the following subsectiori^. So, what we need 
to do is just to expand it with respect to p and perform the Laplace transformation. This 
computation is just a repetition of the matrix model case0 However, as we pointed out in 
the previous section, we should be able to extract more information about the role of the 
resolvent in the gravitational point of view. In the rest of the paper, we will investigate 
this issue. 

3.1 A review of the Drukker-Fiol D3-brane solution 

Let us give a brief review of the D3-brane solution derived in [9], which will be useful for 
later discussions. Their proposal is that the expectation value of the Wilson loop with 
winding number k is given by the amplitude of a D3-brane carrying string charge k and 
attached to the loop on the AdS boundary: 



Here S'ds consists of the Dirac-Born-Infeld action S'dbi and the Wess-Zumino term S'wz 
for a D3-brane: 



^ More precisely, only the leading term J-q, i.e. the leading term of JF_, has been reproduced from the 
D3-branc solution. Note that in the previous section, the pole has been derived by using the explicit form 
of J^o and also the fact that .F_ can be expanded as = J-o(k) + {\/N)!Fi{k) + • • • , but the explicit 
form of higher terms, including Ti are not needed. 

Here we should mention that this argument does not fix the residue of the pole. In order to do this, 
it is necessary to fix a normalization factor in AdS/CFT correspondence up to the subleading order in 
There may exist another difhculty in determining a relative normalization between the cut and the 
pole. This is because, based on the viewpoint discussed at the end of the previous section, these two 
contributions correspond to different setups in the gravity side, namely different perturbative vacua in 
string theory. 




(3.1) 





(3.3) 



(3.4) 
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Here 'P[Ci\ is the pullback of the four- form potential C4 and the D3-brane tension is 
given by T^a = ^27r)3;4g • "^^^ path integral on the right hand side of fl3.ip is taken 
over the D3-brane configurations satisfying appropriate boundary conditions which are 
specified in terms of the Wilson loop. The boundary term S'b is introduced in order to 
flip the boundary conditions for the world volume gauge field and the scalar fields. In the 
semi-classical regime, the path integral may be well estimated by a saddle point value. 

The authors of the paper started with the following metric for the AdSs part in the 
Poincare coordinate: 

r2 

ds^ = — (dy^ + drl + rldijj'^ + drj + rldcj)'^) , (3.5) 

y2 

with the curvature radius L = {AngsNY^^ls = X^^^h, and considered a circle ri = R, 
r2 = on the AdS boundary y = and then discussed the D3-brane solution which is 
attached to the circle. 

The actual D3-brane solution and the evaluation of the action were discussed by using 
the following metric: 

ds'^ = (dr]'^ + cos^ r/rf^/'^ + dp'^ + smh^ p(de^ + sin^ Odcj)'^)) . (3.6) 

sin rjy / 

Here, rj, p and 6 are related to ri, r2 and y by the coordinate transformation 

R cos f] R sinh p sin 6 R sin rj 

cosh p — sinh p cos 9 ' cosh p — sinh p cos 9 ' cosh p — sinh p cos 9 

The range of these new variables are taken to be < p < 00, < 9 < n and < 1] < n /2 . 
They chose tp, p, 9 and (p as the world volume coordinates by = ip^a"^ = p,a^ = 9 and 
= (f) and assumed an symmetric and also ^/'-translational invariant ansatz, rj = r]{p) 
and F^p = F^p{p) . Then Sdbi and Swz are written down as 



f , ,^sin6'sinh p / ^ , 0. , N^sin r? ^„ . 
^DBi = 2N j dpd9^-^^cos^r]{l + r]'') + (2na'y^F^^, (3. 

f cos ?7 sin 6^ sinh^ p / , sinh p — cosh p cos ^ 

i^wz = -2N / dpd9 cos?7 + ?7 sin 77 — — 

J sm rj \ cosh p — smh p cos 9 

and the solution for the equations of motion is given by 



ik\ ( k^\\ 



sin, = K-'si„hp, F,, = ^-^^-^, (k=— j. (3.10) 



Here k is the string charge defined ai 



k = Vil = i I d9d^——, (3.11) 



dF, 



-"^^In this article we follow the notation of [9] for the D3-brane solution except this 11 which differs in 
the factor i. 
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where L is defined as in (13.21) . The solution is attached to the circle on the AdS boundary 
at p = . 

The boundary terms are introduced in order to take account of the appropriate bound- 
ary conditions; the D3-brane carries a fixed number of the string charge /c, and also the 
scalar fields transverse to the AdS horizon should satisfy the Neumann boundary con- 
dition. The latter boundary condition is motivated by the T-duality argument for the 
case with the fundamental string world sheet description instead of the D3-brane [2] . 
The boundary condition for the gauge field is taken into account by adding the following 
boundary term: 

■ (3.12) 



Sa= dtpdOdcj) 

J d{dpA^) 

On the other hand, the relevant boundary term for the scalar fields is the one for the 
radial direction of AdS. Here we need to take care of the fact that the boundary terms for 
the radial coordinate differ for different choices of the coordinates system. In [20] , it was 
argued that a natural choice is to set the Neumann boundary condition for the coordinate 
u = 1/y: 

f dC 

Su= di/jd6d(j) u . (3.13) 

One of the reasons for this choice was that the conjugate momentum of u is finite at the 
AdS boundary, while that of y is infinite. At the same time, this choice is necessary for 
the correct cancellation of divergence. In the following subsections we will see this choice 
is natural in our case as well. 

Summing up all the terms (13.81) . (13. 9p . (13.121) and (13.131) . the classical value of 5*03 + 5'b 
is evaluated as 

+ Sb = ~2N{kVI + k-^ + arcsinh/t) , (3.14) 

which indeed agrees with the matrix model result (12. 5p [9]. Therefore, as mentioned at 
the beginning of this section, by replacing k with k = (k + p)\^/{AN) and performing 
the Laplace transformation for e"*^'^^^'^'^'"'' with respect to p {p <^ N) , the pole at z = 
= a/A^I + can be reproduced. 

In the above discussion, we used the explicit value of the total action with the solution 
substituted, which results in a function of k. As a result, the origin of the pole is not 
very clear. Since the string charge carried by the solution should have been set by the 
boundary condition before solving the equations of motion, we should first change the 
boundary condition as k ^ k + p and then observe a response of the action. It would 
allow us to discuss the Laplace transformation without using the explicit form of the 
solution from the beginning. 

Before doing this analysis, here it is useful to take much care of the boundary condition 
in the above analysis. In the above discussion, we imposed the boundary condition for 
the conjugate momentum of the gauge field as ill = A; . On the other hand, we did not 
mention the condition on the boundary value of the conjugate momentum of u; it is just 
given by inserting the solution. In principle, the boundary condition for the conjugate 
momentum of u should be specified independently before we discuss the solutions or the 
equations of motion. Moreover, based on the spirit of the AdS/CFT correspondence. 
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the boundary condition should be specified in terms of the information contained in the 
Wilson loop. 

In the case of the correspondence between the Wilson loops in fundamental represen- 
tation and string world sheet, a detailed discussion on the boundary condition was given 
in [2]. In the next subsection we consider a boundary condition for the D3-brane. This 
analysis plays an important role in subsection 13. 3[ 



3.2 D3-brane boundary condition imposed by Wilson loop 

As we have mentioned, the boundary condition for the D3-brane should be specified in 
terms of the information carried by the Wilson loop. Before discussing the connection 
with the Wilson loop, let us take a look at what boundary condition is satisfied by the 
actual D3-brane solution fl3.1Up . Near the AdS boundary, the radial coordinate u can be 
approximated by m ~ l/{Rr]). Hence the conjugate momentum of u, at the D3-brane 
boundary, is given by 



d{dpU 

Here we introduce a new coordinate U as 



(3.15) 



p=0 



2na' An 



L2 L2 1 ^ ^ 

U = u = . 3.16 

27ra' 2na'y ^ ' 

The coordinate ?7 is a natural choice for the scalar field, since Li^ /y corresponds to the ra- 
dial coordinate of the fiat six dimensions in the asymptotic region of the original extremal 
black 3-brane solution. 

We define Pjj as the conjugate momentum of U integrated over S^: 

then the boundary condition is simply given by 

li = -ik, Pu = -Rk. (3.18) 

Note that the factor sin6'/(47r) in (13.151) originates from the symmetric property of the 
solution. Hence if we assume symmetry, we lose no information by integrating over the 

In the case of string, the connection between the boundary value of the momentum 
and the information on the Wilson loop is proposed and its validity is supported by an 
argument using the Hamiltonian constraint in [2]. In our present case for the D3-brane, 
we consider the other way round: we examine what kind of conditions such a constraint 
implies on the boundary value of the momenta for the world volume fields, from which 
we can guess appropriate boundary conditions. 

For this purpose, first we note that the D3-brane solution in the previous subsection 
can be locally approximated by the solution corresponding to the straight Wilson line. 
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For later reference, we clarify this point. Let us consider a small region around a point 
on the loop, say = ipo and ri = R, and expand the coordinate aa^n 



= ipo + et , ri = R{1 + ep cos 9) , r2 = Repsm9, y = Ret] . (3.19) 

Here e is a small parameter and we assume that the coordinates t , p , and rj take values 
of order 1 . The angular coordinate 9 can take value in the original range as < < vr . 
Then the AdSs x metric is reduced to the following one: 

ds^ = -^{dif + dt^ + dp^ + p^{d9'^ + sin^ 9d(f)) + L^d^l, (3.20) 

and the Wilson loop is extended along the t-direction. The solution f l3.10p becomes the 
one corresponding to the straight Wilson line in e — > limit. In the following we discuss 
a kind of Hamiltonian constraint only for the case of the straight line. However, from this 
observation, we expect that the analysis can also be applied for a local patch of a Wilson 
loop of an arbitrary (smooth) shape. Here we note that as we will see, when we interpret 
our analysis as the local version of the analysis for a circle, we must take the scale factor 
R and e into account, which is lost in the case of the straight line. 

Moreover, it should be stressed that in the following argument, although we use the 
symmetric ansatz, we do not use information on the classical solution reviewed in 
the previous subsection (nor its reduction to that for the straight line) in the gravity 
side. In our semi-classical analysis, it seems natural to allow any field fluctuation around 
a classical configuration on the right hand side in (13.11) as long as it is consistent with 
the boundary condition. From such a point of view, the symmetric ansatz is just an 
assumption for allowed fiuctuationslll Of course requiring global symmetries including the 
supersymmetries is sometimes so strong that an allowed configuration will be uniquely a 
classical solution, especially when the Wilson loop is highly symmetric as for the straight 
line and the circular loop; see appendix [Bl However here we emphasize that the 
symmetric ansatz still allows a wide range of fluctuations, and the following argument 
should hold for such field fluctuations. 

We use the following coordinates for the AdSsxS^ geometry: 

ds^= (^^y U'{dt' + dp' + p\d9^ + smHd<p^)) + L'^^^ (3.21) 

= Gttdt' + Gf,pdp^ + G0e{d9'' + sin^ ^rf0^) + Guw^ {dU' f (3.22) 
= G xM xi^dX^^ dX^ . (3.23) 

The four directions {t, p, 9, 0) are identified as the four- dimensional space where the gauge 
theory lives. The Wilson loop is located at p = and extended into t-direction. We are 

^^This corresponds to p ^ ep, -q ^ erj. 

It may also be possible that, as a definition of the correspondence, the functional integral on the 
right hand side of p.ip should be constrained in a domain where the world volume fields preserve some 
symmetries of the Wilson loop. In such a case, the symmetric ansatz becomes more important than in 
the case where it is just an assumption. Although it is an important problem to give a precise definition 
for the integration region on the right hand side of (|3.ip . it is beyond the scope of the present article. 
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now discussing the D3-brane configuration which is pinched to the hne p = on the AdS 
boundary, U = oo . We introduce the world volume coordinate a" (a = 1 ~ 4) and, as 
mentioned above, assume the symmetric ansatJ^: 



a 



3 J. _ ^4 



, (3.24) 

t = t{a\a^), p = p{a\a^), U' = U\a\a^), Fu = Fu{c7\ a^) . (3.25) 

We take the coordinate o"^ in such a way that the D3-brane boundary is at cr^ = 0, i.e., 
p(cr^, cr^ = 0) = and U{a^, cr^ = 0) = oo. Then on the boundary = 0, the parameter 

should parameterize the loop on which the D3-brane is attached. 

Now the D3-brane action is given by 

<S'dbi + 'S'wz = J da^da'^Lr,3 = J da^da'^{LDBi + Lwz) , (3.26) 

Ldbi = ^^TMe^J 911922 - gh + (27ra')'^i'2 , (3-27) 
Lwz = -^7iTj,sGeey^G^p{d^td2P - ^ip^st) . (3.28) 

Here the Lagrangian Ld3 and Ldbi.wz are defined with the part integrated. 

We define the conjugate momentum Pxm and 11" for and Aa, respectively as 

dLY)3 a 9Ly)3 , . 

Then, the following identity holds: 

G'\Pt-Vtf + G'''^{Pp-Vpf + G'''''\Pu.f + -^^^{Ii'f9ii = {^nTMe?9ii, (3-30) 

where gu = Gtt{ditY + Gpp{dipY + Guiui^diW^Y . The "momenta" Vxm are defined as 
VxM = (9L-wz/<9(52X*^) , whose explicit forms are given by 

Vt = ATiTj.^Gee^G^pd^p, Vp = -AnTMeVG^pdit . (3.31) 

These do not include any derivative with respect to and thus fl3.30p can be regarded as 
a constraint in the phase space. By using the explicit form of the metric, the constraint 
f l3.30p can be rewritten as 

1 fivr 

= (Pu.-f + (n')2((ai«)2 + (a,p)2) - —N(p,dip-PAt)p' 

- (64^X-'N^(Upr - j^(ny) + ,3.32) 



^''There are several in this AdSsX geometry to be identified with the world volume S^, for 
example, some C S^. Here we consider configurations which are trivial in the direction. This is 
also an assumption. 
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Let us assume that the D3-brane configuration asymptotically satisfies the following 
conditions in the limit cr^ — 0: 

f/^oo, f/p<oo, Pt, P, < Pu,U'ir^ N). (3.33) 

The first condition has already been imposed as a boundary condition, which requires that 
the D3-brane should be attached to the AdS boundary. The second means that the radius 
of the part of the D3-brane is finite at the boundary. The last one refers to the behavior 
of the fields at the boundary with respect to A^. Although the last two inequalities are 
new assumptions, we expect that fl3.33p picks up a reasonable set of fiuctuations. 

Under the conditions (13.331) . the equation (I3.30p . in the limit — 0, is reduced to a 
simple form: 

(Pt;0' + (n')'(5ii)' = 0. (3.34) 

Thus we find that the world volume fields on the D-brane attached to the Wilson loop at 
the AdS boundary are subject to this constraint, at least under the symmetric ansatz. 
Now it is natural to impose the boundary condition on 11^ as 

U\(j\ = 0) = -ik, for each a\ (3.35) 

since the boundary condition imposed by the straight Wilson line, with constant ^* , should 
not depend on the parameter . Then since the conjugate momentum of W and that of 
U , which are defined by (13.291) . are related through P[/» = {U'^/U)Pu on the boundary, 
the equation (13.341) tells us that at cr^ = 0, 

= -ik , Pu^ = Pu&i = -k\dit\ei , Qi = UyU . (3.36) 

Note that there should be a sign ambiguity for Pu. Here we took the minus sign so that 
it will be consistent with the classical solution. This is the only input from the classical 
solution. 

As an application, in order to reinterpret (13.361) in the case of the circular loop, we need 
to take account of the scale factor R and e as pointed out in the discussion after (13.201) . 
Since the first equation of (I3.36P dose not have any uncontracted target indices, it does 
not need to be changed. On the other hand, the second equation should be reinterpreted 
by taking account of the scale factor as [/(straight) = Re x [/(circle) as imphed by the 
last equation in (13.191) . and t = e^^ip . Hence we have 

Pu = —k\dit\ (for straight line) Pu = — /c-R|i9iV'| (for circular loop). (3.37) 

This is precisely the relation we found at (I3.18P (note that the (I3.37P dose not depend on 
the choice of the parameter a^) . 

So far we have considered the local straight line (I3.36P or the circular loop (I3.37p . In 
more general case, we thus expect following boundary conditions: 

H^ = -ik , Pu^ = -k\diX\ei . (3.38) 

Here X^^ is the four-dimensional Cartesian coordinate, i.e., ds\^^^ = Y~'^{dY'^ + dX'^dX^) 
which is identified with the space where the gauge theory lives. 
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In the spirit of the AdS / CFT correspondence, the boundary condition should be given 
in terms of the Wilson loop: 



WkiC) =tTPexp(^j dsf^iA^x" + <!>i\x\9'^^ . (3.39) 

Thus we find that natural boundary conditions are 

X'^{a' =a\s),a^ = 0) = x>'{s), (3.40) 
U\a\(r^ = 0) = -tk, Pu,{a^ = a^{s), = 0) = -k\x{s)\ei{s) , (3.41) 

where we omitted the trivial S^-dependence in (13.401) and cr^(s) determines a relation 
between two parameters of the loop (or the boundary of the world volume) representing 
the same point. Hereafter we assume o"^ = s by using the reparametrization invariance of 
the Wilson loop. We find that the Wilson loop gives the Dirichlet boundary conditions 
fl3.40p on the scalar fields along the world volume direction of the N D3-branes and the 
Neumann boundary conditions (I3.4ip on the gauge field and the scalar fields perpendicular 
to the world volume direction, as expected [9]. 

We would make some remarks. First we notice that in (13.411) the integration over 
in AdSs is done in 11^ and Pfji, which is seen from the definitions (13.261) and (13.291) . 
We expect that this integration is trivial even for loops of arbitrary shape if the loop is 
locally approximated by a straight line with the symmetric ansatz. Next, we note that 
our boundary conditions are natural generalization of the ones given in [2] in the sense 
that (13.401) and the second boundary condition in (13.411) amounts to them in the case 
of the string world sheet, namely when we neglect the a^, cr^-dependence and use the 
Nambu-Goto action. 

Finally let us make a remark on the most important aspect of our boundary conditions. 
By use of the embedding coordinates X^, we can convert the world volume indices of 11" 
into the space-time indices along the world volume like 11'^ = daX^U"' . Since 11° has the 
only non- vanishing component for a = 1, 11^ satisfies the following boundary condition 

= 9iX^n^ = -ikx^" , at = 0, (3.42) 

where we have used (I3.40p . Using the second boundary condition in (I3.4ip . this gives the 
following relation 

(n^'f + {Pu^y = . (3.43) 

It is worth noting that the boundary condition (I3.4ip imposed by the Wilson loop (13.390 
thus corresponds to the BPS condition in [21], i.e., force balance between the electric 
charge 11^ and the deformation of the D3-brane which is characterized by Pu , in the 
case of the spike solution in the flat space. The force balance equation becomes simple 
and symmetric form for our choice of the coordinate, i.e., and W . This may be 
understood by the fact that the coordinate system {X^,X''~^* = 27ia'W} corresponds to 
the Cartesian coordinates of the fiat ten dimensional space in the asymptotic region of 
the black 3-brane solution. Thus it is natural that the equation (13.430 also implies a local 
BPS condition for the Wilson loop. In fact, in gauge theory side, we can introduce a 
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loop of an arbitrary shape in the internal space by replacing \x\d^ in (13.391) with some 
function ?/*(s). For a corresponding D3-brane solution our boundary condition becomes 
the following symmetric form: 

W{a\ = 0) = -ikx^'ia^), PuA<^\ = 0) = -kyi{cT^) , (3.44) 

and then the force balance condition is equivalent to the relation 

x'' = y\ (3.45) 

which is nothing but the local BPS condition for Wilson loops in the gauge theory side 
[2,22]. Therefore we find our boundary condition quite natural because once it is assumed, 
the local BPS conditions in both sides become equivalent!^. 

Another interesting point is that the Gauss law constraint implies that the coordi- 
nates and cx^, i.e. normal direction and tangential direction of the loop, are mutually 
independent. Indeed, from (I3.4ip and (13.441) 



= = g^,i-^kxn, -^kx^ = = = -^^Hk), (3.46) 

and hence da'^/da^ = 0. We would also like to note that (13.441) makes the boundary terms 
(I3.12P and (I3.13P manifestly invariant under the reparametrization0 

3.3 Laplace transformation and matrix model resolvent revis- 
ited 

Let us reconsider the Laplace transformation of the circular Wilson loop taking account 
of the following boundary condition 

n = -ik, Pu = -Rk . (3.47) 

We now evaluate the right hand side of (13. ip under this boundary condition. In a region 
like A ^ 1 the saddle point approximation will be valid, and then we can estimate 
the summation S'ds + S^, by its saddle point value. Since adding the boundary terms 
amounts to the Legendre transformation, this saddle point value is given by a function of 
the boundary values of momenta as (S'ds + Sk) |saddie point = "^(H, Pu) , where 11 and Pu 
express the boundary values of the momenta^ Then the AdS/CFT correspondence for 
the Wilson loop in a semi-classical limit claims 

}_^^^kM\^-ni^^k,-Rk)_ (3.48) 



""^^However this equivalence may be a peculiar feature of AdS^ x background [23] . 

^^It might be interesting to notice that the boundary term p.l2p and p.l3p takes the form of the 
exponent of the U{1) Wilson loop under our boundary conditions, though we should not be confused 
about the gauge field in the super Yang-Mills theory and that on the D-brane world volume. 

Here we still assume the symmetric ansatz. Hence only the conjugate momentum of ri{p) and 
Fp^{p) are considered. 
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As mentioned at the beginning of this section, in order to derive the matrix model resolvent 
associated with the isolated eigenvalue, we change the boundary condition as k —>■ k + p 
and make the Laplace transformation with respect to p , with taking account of the relation 
p -C . Now let us follow simply this procedure without using other information like 
the explicit form of the saddle point configuration in contrast to the derivation described 
below (13.141) . In fact, this approach enables us to clarify a gravitational interpretation of 
the isolated eigenvalue as we will see shortly. 

Recalling that the variation of TZ with respect to change of the boundary value origi- 
nates only from that of the boundary terms [9] aJ^ 



57^(^,Pc;) = [ #5(n)A^ + [ dtps{Pu)u 

J p=o J 



(3.49) 

p=0 



By neglecting the terms which vanish in the limit p <^ N, TZ{—i{k + p), —R{k + p)) is 
rewritten as 



ni-ik, -Rk) -v(^j d'^iiA^ + RU)^ . 



(3.50) 



Then by performing the Laplace transformation, we can derive the following pole term: 

„-n{-ik-Rk) 

,pe-.^e-.-'.'-.^-.--. ^ ^_;,^(,^^,^^^^ . (3.51) 

As a check, we shall plug the solution (I3.10p into this, and then the location of the pole 
reproduces the position of the isolated eigenvalue correctly: 







j dtp{iA^ + RU) = VXVI + ^2 . 



(3.52) 



As we have noted a couple of times in this paper, it is not so striking a result in itself 
that the position of the isolated eigenvalue is reproduced. Here we claim that the position 
of the pole, namely the position of the eigenvalue, has been identified with the electric 
flux Fp^ integrated with respect to p and ip directions of the D3-brane (plus contribution 
from scalar field U) as far as the isolated eigenvalue is concerned. This is exactly what 
we anticipated in the last of section [21 because this value is conjugate to the electric flux, 
namely the string charge, in the sense of fl3.49p and (13.501) . 

Here it is important to recognize that in deriving the gravitational interpretation of 
the eigenvalue (13.521) . our boundary condition (I3.47P plays a crucial role. This is because 
each boundary term (13.120 and (I3.13P diverges at the boundary, but after an appropriate 
regularization, dependence on the regularization cancels out between them as shown in [9] . 
One can easily find that under our boundary condition (I3.47P , the divergence again cancels 
because of the fact that both of boundary values of H and Pu are given in terms of the 
same k. In particular, if we impose the boundary condition only on H, this is not the 
case. Finally, it is interesting to see that from (13.520 the gravitational counterpart of 
the isolated eigenvalue takes the form of the exponent of the U{1) Wilson loop (see the 
footnote at the end of subsection 13. 2p . 



-•^^It is important to notice that this is true because of the saddle point equation. Hence this is not 
generally true for general configurations. 
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4 Conclusions and discussions 



In this paper, we have analyzed the gauge theory Wilson loop winding around a circular 
loop k times, by using a D3-brane carrying k units of string charge in the context of 
AdS/CFT correspondence. It is known that the calculation of the expectation value of 
this Wilson loop, thanks to its symmetry, boils down to considering a Gaussian matrix 
model with an exponential operator insertion. We have then aimed our goal at establishing 
a gravitational interpretation of the eigenvalue in this matrix model based on the D3-brane 
version of the correspondence (11. II) . This point presents a contrast to the preceding papers 
on "bubbling Wilson loops" [12-14]. 

We started with analyzing the resolvent of the Gaussian matrix model using the or- 
thogonal polynomials. In the calculation, we first derive the resolvent in its inverse Laplace 
transformed form, by use of e^^ . Then after the Laplace transformation we obtained the 
usual resolvent. In the inverse Laplace transformed form, we derived the expression which 
is valid to all orders in 1/A^, and we identified a term that corresponds to the eigenvalue 
which is isolated due to a large k ~ 0{N) effect. The remaining terms are responsible 
for the rest of the resolvent. The part corresponding to the isolated eigenvalue has the 
structure in which p is included in the shift of k as tr — >■ tr qC'+p)^^ . In accordance 
with the D3-brane description of the matrix model operator j^tie'''^ , we were naturally 
led to consider the D3-brane with k+p string charge and the Laplace transformation of its 
amplitude. Eventually we identify the position of the isolated eigenvalue, observed as an 
isolated pole of the resolvent, with an integrated flux (and scalar fluctuation) on the D3- 
brane. We therefore succeed in providing, at least in part, a gravitational interpretation 
of the eigenvalue in the Gaussian matrix model. 

As a by-product, we have proposed natural boundary conditions for the D3-brane 
configuration with fluxes in terms of the Wilson loop. These boundary conditions also 
provide a direct relationship between local BPS conditions in the gauge theory side and 
that of the effective theory on the probe D3-brane. 

Let us now discuss the cut of the resolvent. In the present paper, we have mainly 
studied the gravitational description of the isolated pole of the resolvent. It is surely nice 
if the cut of the resolvent can also be discussed based on our gravitational description. 
In the gauge theory side we saw in subsection 12.21 that the cut of the leading semi-circle 
originates from the leading term in w{k,p). In the gravity side, we argued that this 
term corresponds to the configuration with a D3-brane carrying string charge p (-C A^) 
in addition to the one with the string charge k (~ A^) which is now treated as a part of 
the background. So, we may expect that by performing the Laplace transformation of 
the contribution from such configuration, the cut would be reproduced. However, here we 
should recall that the saddle point value of the gravity side gives only the leading term of 
the gauge theory observables. In fact, the saddle point value of the D3-brane action with 
string charge p (^ A^) naively gives the result e^^ . This is just the leading contribution 
of the Bessel function (12.281) in the large 't Hooft coupling limit, and it dose not lead to 
the cut after the Laplace transformation. Hence, it is clear that in order to reproduce the 
cut, we need to take account of the quantum a' correction of the D3-brane with string 
charge p. 
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Another point we should also mention is the case of the correspondence between 
anti-symmetric Wilson loop and D5-brane [10,24]. It would be interesting to examine 
whether our gravitational interpretation is valid also in this case. However the situation 
is not straightforward. This is because the multiplication of the k-th anti-symmetric 
trace operator Tr^j.e*^ and the probe operator tr e^^ dose not include the {k+p)-th anti- 
symmetric trace operator. It seems to suggest that we need to consider the configuration 
with a D3-brane carrying p flux around the background D5-brane with k flux, instead of 
only considering a single D5-brane with {k + p)-fiux. 

Here we also point out a subtlety concerning the relation between a Wilson loop with 
winding number k and that in the k-th symmetric representation. Let Ui be i-th eigenvalue 
of e*^ where M is the matrix variable used in the matrix model analysis in section 12. 2[ 
The Wilson loop with winding number k can be calculated by tr e'^^ = while the 

Wilson loop in the k-th symmetric representation corresponds to 



where ■ ■ ■ includes the other combinations of the powers of the eigenvalues. See [10, 19] 
for the details. Here we note that the decomposition that appeared in our matrix model 
calculation (12.351) is equivalent to the terms here when we write them as w{{k — 1) + 1) + 
w{k — 1, 1) with k being k — 1 and p being 1. It has been argued that the expectation 
values of these two Wilson loops coincide under the large 't Hooft coupling limit, namely 
in (14. ip . compared to the first term, the rests are exponentially suppressed in large A. 
However our explicit calculation shows that w{k) is of order 1/A^ compared to w{k — 1,1). 
Note that k and p have been assumed to be of 0{N) and 0{1) respectively, and then the 
calculation is still valid. The equivalence thus holds only when the large-A limit overcomes 
the difference in 1/A^, that is, A needs to be larger than log A^. This is therefore out of the 
usual limit such as first taking N —>■ oo with fixed A and then taking A to be large. Note 
that in this article we employ only the fact that in the strong coupling limit the D3-brane 
solution with fc-fiux agrees with tre'^*'^, and then this subtlety does not matter. 

Finally, it would be also important and interesting future work to clarify the relation 
between our viewpoint and the bubbling picture [13,14]. We believe our viewpoint leads 
to a deep understanding of the gravitational interpretation of the eigenvalue, or moreover 
the connection between gauge theory and gravity. 
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A Two loop correlator in Gaussian matrix model 

In this appendix we compute the "two loop correlator" of the Wilson loop: 

W{v). = yJ c^Mltre-^^ltre'^^^e-^-^^ , (A.l) 



1 



Z^^ = j dM-tie^^'e-'-^''^'\ (A.2) 

in order to examine the resolvent in the presence of the Wilson loop with a large winding 
number. For this purpose, we have to evaluate (lA.ip for finite N as noticed below (12.331) . 



Therefore, we calculate (lA.ip by means of the orthogonal polynomials. We first change 
the variables as 

^'^l/^*. P'=\/^P. ^''=y?^'. 

then W{p)k becomes 

W{p), = i I rfM'ltre^'^^'ltre^'^'e-^-^" , (A.4) 

Z, = ' Z, = l dM'^tre'^'^^'e-'^'^'^" . (A.5) 

Using the Hermite polynomial Pi{m) satisfying 

dme-^"'^Pi{m)Pj{m) = hi5ij, hi = V2^i\ , Pi{m) = m' + 0{m'~^) , (A.6) 



the Vandermonde determinant can be written as 

, 2 

tP,-^i(mi) ' 



A{mf = (detPj^i{mi)^ = ^ sgn((Tr) JJ P^(^)_i(m^)P^(^)_i(m^), (A. 7) 

^ 'or t 

and the relevant integral is given by 
ZkW{v\ = I \[dm^e-"^^-Y.^gn^''^)\[P„(i)_,{mi)P,^i)_,^^^^ 



'^3 

(A. 



We decompose the sum over j in (lA.Sp into two types according to i = j oi i ^ j 

(lAl) 

+ E / n^"^^^"'"^ n^-(0-i(^')^.(0-iK)e(^'+''^"^'. (A.IO) 



i k 
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Each term flA.91) and flA.lOl) corresponds to Zkw{k,p) and Zkw{k + p), respectively where 
w{k,p) and w{k + p) are given in fl2.35l) . 
The first term can be rewritten as 

(MD 

= Yl J dmke~^'^^*^Pa{k)-i{'mk)Pr{k)-i{mk) 

x/.^,e-^™?P.,,_.K.)P.OV.K)e^'- 

X J dm,,e-^-?P^(,)_i(m,)P.(,)_iK)e^'''"' (A.ll) 
^ ]^ XI ^S^*^^^) 5Z 5Z ( n ^^{kMk)haik)-iyaU)-iMj)-dk')Lii)-iMi)-iiP^^ ■ (A- 12) 

err 

Here we defined 



^km 



(A.13) 



On the other hand the second term flA.lOp can be rewritten as 

err i j^i 

X j rfm,e-^'"?P.(,)_iK)P.(,)_i(m,)e(P'+'=')"^> 

^ ^sgn((Tr) XI (n^'^0>(j)^'^0)-i) J dmie-^"^^ P„^,yi{mi)Pr(iyi{mi)e^P'^''''>"'* 

(A.15) 
(A. 16) 



(A. 14) 



E E ( n ^'^a)-i)^<^«-i(p' + ^'), 



cr j j^i 



with /j(/i;) = Ii^i{k) . From (lA.lSP to (]A.16p . we used the fact that the Kronecker delta 
Wj^i ^a-{j)T{j) implies that two permutations a and r are identical. 

Next we perform the integral in (lA.lSp . We first note that from the generating function 
of the Hermite polynomial 



^ — ^ 2! 



we have 



i=0 



t=o 



(A.17) 
(A. 18) 



Plugging this equation into flA.lSp . we have 



Iij{k) = I rfme-^™'+'="(9^e*"-^9^e'™-^ 



s=t=Q 
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dldle 2 6 2 



s=t=0 



=t=o 



(A.19) 



t=o 



By differentiating the generating function of Laguerre polynomial L^-"^ {x 



(l+t)"e-* = ^Lf-')(x)f 

j=0 



(A.20) 



with respect to t we obtain 



dt 



-{l + tre 



a^—xt 



(A.21) 



By using this equation with replacing x —k^ and t t/k, we can rewrite (]A.19|) as 

Ii,,{k) = hie^'''y-'L^r'\-k^) ■ (A.22) 

Notice that by definition 

h^j{k) = Ij,i{k) , (A.23) 
which can also be proved explicitly from flA.221) . and 



Substituting (1X221) into ([02D yields 

(jAU) 

^ ^ sgn(ar) ( J] ^-(k)r{k)K(k)-i) 



LAx) = LfHx). 



(A.24) 



^sgn(ar)(^ H ^-ik)r(k)) 

ar I i jj^i k^i,j 

-k'^)P 

I ar i j^i k^i,j 

^ J^fr{j)-a{j)^{T{j)-a{j))^ _ j,,2^^,r(i)-<T(j) ^ (r(i)-<T{i)) ^ _ ^,2 



L 



^,.u>-o,j,,^_j^,.^^^^^^^^ (A.25) 
Here the Kronecker delta in the last expression implies that only two cases are possible: 
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1. cr = r, 

2. a{k) = T{k) for k ^ i,j and a{i) = r{j), a{j) = T{i). 
According to this, we get 

z 

X (L,-,(-«:'')i.-,(-p'') - *"-'4i7'(-«=")p""'iSir'(-p'')) 

(n"-) 

X - k"-'Lf:l\-k'V-^L^{-j/')) , (A.26) 

where in the last step we have used an identity of the Laguerre polynomial, Li^~^^\x) = 
On the other hand, using (lA.24p in (lA.16p . we obtain 

a i j 
j a i 

=|5<=''''""''(nv.)Ei-.(-(p'+*'« 

= ^e^^^'^^'^Xn + (A.27) 

i 

From flA.261) and flA.271) . flA.81) can be rewritten as 
l.H-'(rt.=§ei<'"+'")(n'.,-0 
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+ e''^'L^^U-{k'+pr)). (A.28) 

Moreover, by setting p = in (lA.SP and noting W{0)k = 1, it is easy to see that is 
nothing but (lA.lOp with p' = multiphed by A^. Thus from (]A.27p 

ei^'^(nVi)^S^^-i(-^'')- (A.29) 



a" 



Therefore we find that 



( IA.30p and (]A.3ip come from (1A.9P and (lA.lOp . so give the exphcit form of w{k,p) and 
w{k + p) in fl2.35p . respectively. Notice that so far we have not made any approximation 
and hence this equation is exact and holds for any finite A^. Therefore we can expect its 
useful application elsewhere. 



A.l The evaluation of the second term in w{k,p) 

In this subsection, we evaluate the large- A^ limit of the second term in ( ]A.30p and derive 
the corresponding term in the resolvent. We consider the numerator of the second term: 

N-l 

A=Y,F,,{k',p'), (A.32) 

i,j=0 

Fij{k',p') = e^k"-^Lf'\-k'^)e^p'^-'L^r'\-P^) (A.33) 
= hj'l^,{k')h;'l,,{p'), (A.34) 

where we have multiplied the additional factor e 2 and, for notational simplicity, we have 
changed z — > z + 1, j ^ j + 1. hi and lij are defined as flA.6P and flA.13p . respectively. 

Let us start with the large-A^ expansion of Ij,i{k') for arbitrary z, j and finite k. 
For this purpose, we recall its integral representation flA.13P and rescale the variable as 
m = ml^~N such that the Gaussian potential becomes proportional to A^; 

.(A;') = A^"^ y dme-^'^'+'^''^Pj{m)Pi{m). (A.35) 

Here k' = \fNk' and we have also rescaled the Hermite polynomials as Pi{ni) (m) 
to make it satisfy the normalization Pi{m) = m* + - ■ ■ . The large- A^ behavior of the generic 
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orthogonal polynomial in this normalization is addressed in [18]. Because of the large lin- 
ear term k'rh = 2KNrh, integral flA.35p has a saddle point m = at a non-oscillating 
region ^y4i/N < for the Hermite polynomial Piim) (see [18] for generic behavior of 
the orthogonal polynomial). Hence we concentrate on the expression for that region: 



P,(m) = exp i^N jj d^\ogk^^\m,0 + -logA;W(m, ^) - -logg(m, ^) + 0{N~' 

(A.36) 



where 



m 



+ _ 4^ 



,0 = ^^2-4^. 



By using this expression, we find the saddle point of the integral flA.35P as 



(A.37) 



(A.38) 



Evaluating ( ]A.35p semi-classically including the Gaussian integral around the saddle point, 
we obtain the following large- behavior 



i^'^.+.log 



— — h J log — 



2 2 



N 



^— - ^ log(fc'm,) + ^ log(27r) + 0{N-') 



(A.39) 



Next we study the large- behavior of Iij{p') for finite p = p' ^/WJX. In this case, 
the linear term p'fh = p'm is small and we need to consider an oscillating region. Here, 
instead of studying it, we derive the large- A^ behavior of Iij{p') by solving the differential 

equation satisfied by f{i,j,p') = e^L. 



p 



(A.40) 



The first term {X/AN)p'^ in the round bracket is small compared to the second term in 
the same bracket. By neglecting this term, we can reduce the the equation to Bessel's 
differential equation and we obtain 



kjiP) = v27rmax(z,j)! 



(A.41) 



where A = X{i + j + 1)/2N. max(z, j) represents the larger value and I\i-j\ is the modified 
Bessel function. The overall constant can be fixed by comparing the asymptotic behavior 
of the Bessel function with that of the Laguerre polynomial in the small p limit. This 
expression is again valid for any values of i and j. 



29 



Now we discuss the summation (1A.32P using these ingredients. We should be careful 



that the subleading 0{N~^) terms in flA.39p and flA.4ip possibly contribute to the leading 
behavior of A, since A is defined by the summation over A^^ terms. However, as we will 
see later, only the terms with N — i = 0{1) and N — j = 0{1) actually contribute to 
the second term of (lA.30p . Then the effective number of terms being summed is of order 
(9(1), and we can consistently neglect the 0{N~^) terms. In the following, we call the 
range N — i = 0{1), N — j = 0{1) "C". We first concentrate on the terms in the range 
C, and then we study the contributions of the terms outside C. 

In order to discuss the terms in C, we change variables as m = 2N — {i + j) and 
d = i — j and expand flA.391) and flA.411) with assuming m = 0{1), d = 0{1). Then we 



obtain the following simple expression for Fi j{k',p') 



2N (^kVk^ + 1 + \og{VK^ + 1 + /«))+ log(VK2 + 1 + k) (A.42) 

'l^l^^^^ (l + 0(Ar-)). (A.43) 



^ loginVn^ + 1) - ^ log(87rAr) 



We set the range of summation as = 0, ±1, ±2, ■ ■ ■ , ±D and m = \d\ + 2, \d\ + 
4, ■ ■ ■ ,\d\ + 2M for a fixed value of d. Here M and D define the upper limits of the 
summation which should be large but still of 0{N^) in the present assumption of the 
range C. In the large- limit, both of the upper limits can be consistently taken to be 
infinity. After the m-summation, we obtain the following large- A^ expression for A: 

D M 

A=J2 E F,Ak',p') + {■■■) (A.44) 

d=-D {m-\d\)/2=l 



Here for simplicity we have taken the limit D oo, M ^ oo which is consistent in the 
large- A^ limit as noted above. The extracted factor eVL^^^(— A;'^) is exactly canceled by 

the denominator of ( 1A.30I) taking account of the extra factor e". 

The dots (■ ■ ■ ) in (1A.45I) express the terms outside C. We will shortly see that these 
terms actually dose not give finite contribution to the resolvent in the large- A^ limit. 
Before that we shall derive the resolvent corresponding to the first term of flA.45P by 
performing the Laplace transformation (recall the overall factor — 1/A^ in the second term 
of (jAJQl): 

2ndterml ; Nj^ ^ /-^ (.^^l^^^\\d\ ^ ^ 




30 



Therefore, the only singularity of -R2nd term(^) is the cut —\/~X < z < \/~X. In fact, it is 
easy to see that flA.48p is just another form of fl2.18p . By checking the residue of the 
pole at infinity it is clear that the resolvent corresponds to "minus one" eigenvalue, i.e., it 
subtracts one eigenvalue from the leading semi-circle distribution. Hence we have derived 
the resolvent (12.18^ without assuming the conditions mentioned after ( I2.17p . 

Finally let us show that the terms outside C do not contribute to the resolvent. In fact 
these terms are exponentially suppressed with respect to compared to the terms in the 
range C. Since such exponential suppressions will not be compensated even by summing 
up all the N'^ terms, we can conclude that these terms do not contribute in the large- 
limit. 

In order to clarify this point, it is helpful to start with the following simple ordering 
property of Fij{k',p'): 

F,,,(fc',p')<i^m,i+i(A:',p), (A.49) 

which follows directly from the definition of the Laguerre polynomial: x^) = 

Yll=o j^i-r^^^ Here jCi-r is a binomial coefficient and we define jCi-r = for 
j — i + r < 0. Because of this ordering property and also the symmetric property 
Fij{k',p') = Fj^i{k',p'), it is sufficient to show the following two claims: (I) among the 
terms Fj Ar_i (/;;', p'), those in C are exponentially larger than the other terms, and (II) 
among the terms Fi j{k',p') with i — j = 0{1), those in C are exponentially larger than 
the other terms. By showing these two claims and considering the ordering property 
flA.49p . one can easily see that the terms in C are exponentially larger than all the other 
terms. 

Let us begin with the first claim (I). We use the expressions dAJH) and ([QTIl with 
setting j = A^ — 1. Keeping terms of up to 0{N), we have 



^(e- 1) logAr+ ^(1 -0 + ^R + Nhg^ + N^\og^ 



\N-i{k\p') = exp 

+ N{1 - logp' - AT^log^ - N{1 - log (1 - + 0{\ogN) 



(A.50) 



where ^ = i/N and K± and R are defined as 



^ + 2K±^{l-i), R = .JIQk^ + 8^2(1 + ^) + (1 _ ^)2. (A.51) 



2k 2k 

The first and the second line is the asymptotic form of h'^_^I]^_i^i{k') and h~^IiN-i{p'), 
respectively. We have used the asymptotic form of the Bessel function which follows from 
the saddle point approximation of the integral: 

=p'^-^-^^^^^-j^ ' \^^^^.{pVA){l + 0{N-')) (A.52) 
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where u = N — 1 — i. The second hne of ( 1A.50I) follows by assuming z/ ^ 1, and applying 
the saddle point approximation for the integral and also Stirling's formula for r(z/ + 1/2). 
For the range u = 0{1), in the sense of large- limit, both of these approximations break 
down. However, even in this region, we can find the regime where u is still large compared 
to 1 and also p\/X. In such a regime, the asymptotic behavior flA.501) is still valid. In the 
second line of flA.SOp . we have also used Stirling's formula for i\. It is easy to check that 
the following arguments are still intact even if we consider the regime i = 0{1) seriously. 

The dominant ^-dependence of f lA.SOp arises from the first term on the first line and the 
terms on the second line. Let us take the derivative of the summation of these dominant 
terms with respect to ^: 

(9/1 \ A^'^ 

— Ar(-eiogAr-eiogp'-eioge-(l-01og(l-0) = N log ^—j;^ + 0{N) . {AM) 

In the right hand side, we introduced 7 by A^ — i = A^''', i.e., 1 — = A^'''"-'^ and used 
p' ^ A^~^/^. From this expression, we see that for the range < 7 < 1, Fi^^-iik' ,p') 
increases quite rapidly with respect to C, as N"'^^. 

Let us also check the ^-dependence for the range 1 — ^ = 0{N~^), i.e., N — i = 0{1). 
For this purpose we expand flA.SOp with respect to the small parameter e = 1 — ^. Then 
we have 

Fi^N-i{k',p') = exp A^|^-elog(K + Vl + /«2)-ielog A^+elogp'+e-elogeH ^ (A.55) 

Here we have written down only the leading e-dependence. So, dots (■ ■ ■ ) include the 
subleading terms and also the leading but e-independent terms. By taking the derivative 
of the exponent with respect to = 1 — e, we obtain the following result: 

N log(/€ + Vk^ + 1) + N log (^-^Ne^ . (A.56) 

So, for the range e > p\/X/2N, the function Fj^Ar_i(fc',p') increases with ^ at least as 
g5iviog{/t+\//t^+i)_ ^]^g 0(logA^) terms in ( jA.501) do not change this conclusion. Hence we 
have shown the first claim (I). 

Next we turn to the second claim (II), i.e., the parameter range i— j = 0{1). For this 
range, h~^Iij{p') behaves, at most, like power in A^ as can be seen from (]A.4ip . Hence 
the dominant behavior is determined by the following asymptotic form of (]A.39P : 



e^k"-'L^;~'\-k'^) = exp (^2N(^K^/^^^ + C\og '^^ ^J^^^ ^+0{\ogN)^, (A.57) 

where ( = j/N. The dependence on i — j appears only in the subleading terms. By 
taking derivative of the leading exponent with respect to (, we find that Fi j{k',p') in this 
range grows at least as Fij{k',p') ~ gSA^c iog{K+VK^+i)^ This shows that the second claim 
(II) indeed holds. 
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B BPS conditions and the uniqueness of the solution 



In AdS/CFT correspondence, one of the most important guiding principles to find out 
a corresponding object to its fiolographic counterpart is the symmetry preserved by the 
object. For example, the circular and the straight line Wilson loops preserve S'L(2,M) x 
S0{3) symmetry as a part of the Euclidean conformal group 5*0(5,1) of the A/" = 4 
super Yang-Mills (SYM) theory. In accordance with this preserved symmetry, the D3- 
brane solutions found in [9] have the structure of AdS2X whose isometry is indeed 

The circular loop and the straight line also preserve a part of the 
supersymmetry of SYM, namely they are BPS objects, and therefore corresponding D3- 
brane solutions should preserve some of global supersymmetry of type IIB supergravity 
and, as also shown in [9], so do they. 

In this appendix, we see that the BPS condition for D3-brane solutions, together with 
the S^ symmetric ansatz, suffices to determine the solution (virtually) uniquely, at least 
in the circular loop and the straight line cases. In [9], the authors checked that, for the 
straight line, the BPS equation is satisfied by their solution. Here we do not assume the 
D3-brane solution, but just postulate the S^ symmetric ansatz r] = rj{p) and Fp^{p), and 
observe how the BPS condition restricts the form of the solution. 

We examine the circular loop case with a brief summary of the BPS condition for 
D-brane solutions. We start with the Euclidean AdS^ metric (13.61) . 

ds^ = — ^ (dr]^ + cos^ r] dip'^ + dp^ + sinh^ p {dO^ + sin^ 9 d(f)) , (B.l) 
sin rj 

and the S^ part is omitted in the analysis since the solutions we are interested in are 
trivial on S^. By a coordinate transformation, this metric is mapped into 

ds^ = ^ [dy^ + de + dr^ + r^{d9^ + sin^ Odct?)) , (B.2) 

where Wick rotation to Lorentzian metric, tE = it, is easily understood. One can see 
that this Wick rotated Lorentzian metric is mapped into the metric (IB.ip with replaced 
with ipE = ci-nd then we can consider ipE = iip as 'Wick rotation" of our metric. From 
fIB.ip . we define the vielbeins. 



= dr] , = cos'qdtjj , = dp, 

sin r] sin r] sin r] 

= sinh pd6 , = sinh p sin 6d(j) . (B-3) 

sin r] sin i] 



We then consider the Killing spinor equation. 



D, + ^T,^, ) e = , (B.4) 



-"^^The small fluctuation around the loop also obeys the classification with respect to this symmetry. 
See, for example, [25]. 
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where = ri2345 and the gamma matrices in the flat tangent space are denoted by capital 
Ta and those in the curved background by 7^ = e^r^, and e is a complex combination of 
two chiral Majorana-Weyl spinors of type IIB supergravity. We find the Killing spinor of 
this background, 

e =tan^/2 (^|) (M+er + M_e2 ) + i cot^/' (^|) ri2(M+er - M^e^) , (B.5) 

where are constant spinors satisfying Fej^a = ~^i,2j ^ = r2345, and 

M±{ijj, p, e, (f)) = e^^'^e^^^^^ef^^^e^r^ . (b.6) 

The BPS condition for a D-brane solution is given by a compatibility condition of 
the K-symmetry of the D-brane action and a part of the global supersymmetry of the 
background. See, for example, [26]. The K-symmetry projector with Lorentzian signature 
is defined by 

rfP+i^r = -e-^£^i^,e^AX|^^^, (B.7) 

x = er(2„)X"/, (B.8) 

r(n) A • • • A dC"d,,X'^' ■ ■ ■ a,„X'^"7^,..^„ , (B.9) 

77/. 

where ^ denotes the world volume coordinates, X^^{^) are the embedding coordinates, 
— 27ra'F is the two-form field strength of the world volume C/(l) gauge field, and K 
and / act on spinors as Kip — lip — —iip. This projector is traceless and equipotent, 

Trr = 0, r^^i, (B.IO) 

and the D-brane action enjoys the «;-symmetry 

seiO^{i + m)HO, (B.ii) 

where 9 is the fermionic partner of X^^. 9 is fermionic coordinate of the target space-time 
as well and then it transforms under the space-time supersymmetry as S9 — e. Together 
with the K-symmetry, 9 transforms as 

59 ^ (l + r)K + e, (B.12) 

and therefore if a constant e satisfies 

(l-r)e = 0, (B.13) 

then this global space-time supersymmetry is compatible with the K-symmetry, hence the 
embedded D-brane is BPS. 

We choose p, ip, 9, <p as the world volume coordinates as before, and then with the S^ 
symmetric ansatz, the projector takes the following form after Wick rotation. 



cosr/(l + r/Tia) - 2T^a' ^j^F^pT^^K 



n, (B.14) 
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where 



f =ycos2r^(l + r/'2) + (2vra')^-^i^|p. (B.15) 
Now we try to solve the BPS equation 

Te = e. (B.16) 

First we rewrite the Kilhng spinor flB.5p . 



where 



e=T(r/)(M+er+M_e2) 

+ icoshpC(r/)(M+ri2er - M_Ti2e2) 

+ smhpC{'n){M+g{e,cP)T,e^ + M^giO, (|))T^e,) , (B.17) 



cos 6 + T4 sin 6 cos + Fs sin 6 sin , (B.18) 



and C{r]) = cot^''^ (^^ ^T{r]) = tan^/^ (|) have been introduced for simphcity. The action 
of the projector F on the Kilhng spinor is 

if{r],F)re = - i cosr/r(r/)(M+e^ +M_e2) 

— cos ?7 (cosh pC(r7) + 77' sinhpT(77)) {M^T i2e^ — M_T 1262) 

+ i cos r] {sinh pC {t]) + r]' cosh pT(r/)) {M+g{e, (f))Tie^ + M^g{e, (f))T 162) 

+ t( cosr]C{v){M+g{e, 0)F2er - M^g{e, 0)F2e2 ) 

+ ,2W^F^p (-T(r/) (M+(7(0, 0)F2f + M„^?(^, 0)F2f 62-*) 

~tC{7]) cosh p{M^g{e, 0)Fif er* - M^g{e, 0)Fif £3 *) 

+C(r/) sinhp(M+FiF2fer* + M_FiF2f 62*)) . (B.19) 

In order for a solution to exist, €^2 have to be related to their complex conjugates ej"2 in 
a certain way, say, 

Fi,2fer;2CxFi,2er,2, (B.20) 

and at this stage we need to consider all possible combinations of 1, 2 indices. Since 
and €2, and their complex conjugations, always appear with M+ and M_ respectively, 
flipping the index, say like e^* ^ €2 , is not allowed. So there remain two possibilities: 

Case I: Fife]"* = «iF2e^ , Fif * = a2T2^2 where ai,2 £ C 

First by looking at the signature of eg , ai = 02 is concluded, and thus we take a = ai = a2 
and obtain 



z/(F-l)e = 



i cosTjTivi) — ili^cJ — —^F^paCirf) sinh p — ifTirf) 

Lj 
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+ 



+ 



+ 



= 0. 



— cos 77(cosh pC {rj) + rj' sinh pTijf) ) + / cosh pC {rj) 

i cos?7(sinh pCirf) + -q' cosh pT{r])) + i2TTa' — :jrj^F^paT{ri) — if sinh pCirf) 

ij 

X (M+^(^,0)rier +M_(7(0,0)rie2) 

7' COS 77(7(77) + 2Tia' — -^F^paC{r]) coshp 

X {M+g{e,<P)r2e^ -M_g{d,<f))r2e^) 



(B.22) 



(B.23) 



(B.24) 



(B.25) 



Each coefficient of e terms has to vanish independently, and by ehminating F^p and f{p) 
from these four equations we have 



T]' cot T] = coth p 



(B.26) 



which can be integrated to be 

sin 1] = sinh p. (B.27) 

Here is a constant of integration. By inserting this solution for 
follows from the term flB.24p . we have 



the condition which 



-1 



2iTa' sinh^ p ' 



(B.28) 



and also from (lB.22p . we have 



f{p) = Jl + l{l + a-^) = l 



(B.29) 



which leads a = ±i. It is easy to see that these solve all conditions. 
Thus a half BPS solutions in this case are 

sin?7 =K~^ sinh p , 



2 — ' 

27ra' sinh p 



(B.30) 
(B.31) 
(B.32) 



where the signatures are taken to be same. The solution found in [9] corresponds to the 
plus sign. In order to determine k, one needs to solve the equations of motion for the 
D3-brane action, and we have already known that these BPS solutions solve the equations 
of motion as well. 
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Case II: Te^* = aie^ , Teg* = 02^2 where ai^2 ^ C 
The analysis goes in parallel with the case I, and one finds 

/(p) = -cosr/, (B.33) 

which does not have the solution since < rj < it /2. So this projection does not provide 
a BPS solution. 

We therefore conclude that the BPS condition for D3-brane solutions with the S^ 
symmetric ansatz is sufficient to determine the classical solution. In the main part of 
the text, we have implicitly assumed that as for a small deformation of the boundary 
condition k —>■ k + p there exists a unique classical solution associated with the new 
boundary condition. The result of this appendix justifies this prescription, since the 
solution depends only on the parameter n. 
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